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The volume of separable states is 
super-doubly-exponentially small 



o 

O ; Stanislaw Szarek (Paris and Cleveland) 

Abstract 

\ In this note we give sharp estimates on the volume of the set of sep- 

arable states on N qubits. In particular, the magnitude of the "effective 
radius" of that set in the sense of volume is determined up to a factor 
which is a (small) power of N, and thus precisely on the scale of powers 
■ of its dimension. Additionally, one of the appendices contains new sharp 

estimates (by known methods) for the expected values of norms of the 
GUE random matrices. We employ standard tools of classical convexity, 



■ high-dimensional probability and geometry of Banach spaces. 



Let H = H N := (C 2 )® w be the JV-fold tensor power of C 2 and denote by d 



2 its dimension. In this note we investigate the structure of the set T> = T>^ = 
V(Hn) of states on B(TCn) and, in particular, of its subset S = Sn consisting of 
(mixtures of) separable states. We recall that when T> is identified with the set 
of density matrices {p G B(TC) : p is positive semi-definite and trp = 1}, then 



S = conv{pi <g> p 2 ® . . . <g> p N : pj eV(C 2 ),j = 1,2,..., N}. 

Above and in what follows, we skip the subscript N whenever its value is clear 
from the context. We emphasize that separability of a state on B(TC) is not an 
intrinsic property of the Hilbert space TC or the algebra B(H); it does depend on 
the particular decomposition of H as a tensor product of (smaller) Hilbert spaces. 

The question of the size of S and, particularly, of its relative size as a subset 
of V was raised in [21] and further investigated, among others, in |U UH UH E] 
(see also the survey [15J). One of the parameters that have been studied was the 
maximal size of homothetic images of D contained in S. More precisely, one asks 
for which values of e (say, with e > 0) we have 

eV + (1 - e)I d /d C S, (1) 
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where Id stands here for the identity matrix in d dimensions; in the present 
context, Id/d is referred to as "the maximally mixed state." Alternatively, one 
considers inclusions of type (JTJ) with T> replaced by the appropriate Euclidean 
(Hilbert-Schmidt) ball B. The bounds obtained to date show that, in both cases, 
the optimal (largest) value of e is (asymptotically, as d — > oo) of order contained 
between d~ l and d~ 3 ^ 2 . While clarifying the situation somewhat, all these results 
leave open the question of the precise asymptotic order of the "in-radii" of S on 
the power scale in dim<S = dim 22 = d 2 — 1 as well as the issue of the "size" of 
S when measured by global invariants such as volume. In the latter direction we 
obtain here the following bounds 

c fvo\S\ 1/dimS C(logdloglogd) 1 / 2 

where c, C > are universal effectively computable numerical constants and a := 
log 2 (27/16)/8 « 0.094361 (in other words, d l ' 2+a = 3 3N / 8 ). The "effective radius" 
of S in the sense of volume is thus precisely determined on the scale of powers of 
d. Since complexity of a set can often be estimated using volumetric methods (see 
[T3*| for a modern exposition of this circle of ideas), this goes a long way towards 
the ability to compare complexities of S and of P - even though the so-called 
Bures metric and the related volume may be more appropriate measures of size 
in the present context (see J3 )• In what follows we shall present the main line 
of the argument leading to (J2J), relegating to Appendices the discussion of some 
peripheral issues as well as the description of results and concepts from convexity 
and geometry of Banach spaces that are being used. We refer to jH UHl IE] 
for a more professional exposition of the relevance of separability to quantum 
computation in general and to NMR computing in particular. 

For comparison, we note that the results of 4J and [6 j implied lower esti- 
mates on (vol<S/vol£>) 1/dim ' s which were of order d 13 with, respectively, (3 = 
Iogl0/log4 pa 1.660964 and (3 = 1. By contrast, no non-trivial ■upper estimates 
on the volume of S were apparently available, except in very low dimensions. In 
the opposite direction, the expression on the right hand side of (J2J) yields upper 
estimates on the e's that may work in and related inclusions. However, (J2J) 
does not generally improve on the most recent prior results of that type. One 
case when it does is by giving an upper bound on the radius of a Euclidean ball 
that may be contained in S which is tighter (roughly, by a factor of d a ) than 
the rather elementary 0(<i _1 ) estimate; see © and Appendix H for more explicit 
statements in this regard and for more comments. Here we will just mention that 
our method does not give - at least without any additional work - any explicit 
state that constitutes an obstruction to the inclusion eB + (1 — ejl^/d C S for 
e = o(rf~ 1 ), and that our results suggest that it may be more appropriate to com- 



2 



pare S to an ellipsoid which is substantially different from the one induced by the 
Hilbert- Schmidt norm (see the paragraph containing ()lL)jl). 

Since it is conceivable that the inequalities (J2J) may be of interest not just 
asymptotically, but also for some specific "moderately large" values of N, we 
put some effort into obtaining reasonable (but certainly not optimal) values of 
the numerical constants. Our main argument gives c = 1/4 and shows that 
(J2J) holds with 4(iVlog 2 (4N)) 1/2 = 4(log 2 d log 2 (41og 2 d)) 1/2 in the numerator of 
its third member. A slightly more precise (and more tedious) calculation yields 
c = \/ e/8ir pa 0.32887; see the comments following (jllj) and Appendix E. It is 
also easy to follow the argument and to obtain somewhat sharper estimates for 
specific values of N, which may be of interest, e.g., in the context of a threshold 
of 23 mentioned in Such improvements are sketched in Appendix G leading 
to a non-trivial (i.e., < 1) bound on (vol5/volP) 1//dim5 starting with N = 6 
(by contrast, 4(iVlog 2 (4iV)) 1/2 /d 1/2+a < 1 iff N > 8). Likewise, tighter bounds 
can be obtained if one is only interested in large N; for example, one may have 
c = Cn -»• e 3 / 4 /^ « 0.84456 and C = C N -> e 1 / 4 - s /2/log2 w 2.1811, see 
Appendix E. Finally, our methods allow analyzing separable states on tensor 
products involving spaces C k with k > 2, leading to non-trivial but not definitive 
results; some remarks to that effect are presented in Appendix I. 

Instead of working directly with D and S, we shall consider their respective 
symmetrizations 

A := conv(-X? U V), £ := conv(-S U «S), (3) 

where all sets are thought of as being contained in the real <i 2 -dimensional vector 
space of self-adjoint elements of B(TC) (further identifiable with A4 d a , the space 
of d x d complex Hermitian matrices). We do that because, firstly, the geometry 
of symmetric convex sets is much better understood than that of the general 
ones and, secondly, the specific symmetric sets A and £ are familiar objects 
in geometry of Banach spaces, which allows us to refer to known concepts and 
results. In Appendix D we indicate how one can treat directly T> and S without 
passing to symmetrizations; however, this yields only a very small improvement 
in the constants c, C in (j2J at the price of obscuring somewhat the argument. 

It is readily verified that A consists exactly of those (Hermitian) elements 
of 13(H) whose trace class norm is < 1. Equivalently, A is the unit ball of the 
space Cf := (Mf, || ■ ||i), where, for p G [l,oo), \\A\\ P := (ti(A* Af' 2 ) 1 ^ is the 
Schatten-von Neumann p-norm of the matrix A and Ai s d a stands for the space of 
d x d Hermitian matrices. A similar argument shows that E is the unit ball of 
the iVth projective tensor power of C\ (in the sense of the Banach space theory, 
see Appendix B). We shall denote the corresponding norm on A4 s d a by || • For 
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future reference, we note that in the above notation || • corresponds to || • || op , 
the usual norm of a matrix as an operator on the Euclidean space. We also point 
out that while in this note we focus on (IR-linear) spaces of Hermitian matrices 
and operators, the Schatten-von Neumann classes are most often defined in 
the literature to include all (be it real or complex) scalar matrices and not just 
the self-adjoint ones. 

The plan of the rest of the argument is as follows. First, using classical 
general results from convexity, we relate the volumes of A and £ to those of 
T> and S. Next, we obtain two-sided estimates for vol A and volS, which are 
most conveniently described using the following concept: if if is a subset of an 
n-dimensional Euclidean space with the unit ball B, we call (vol K/volB) 1 ^ 
the volume radius of K. [As hinted earlier, in the present context the Euclidean 
structure is determined by the 2-norm defined above, also often called the Hilbert- 
Schmidt norm or the Frobenius norm, and the inner product is (u,v) = truv.} 
Equivalently, the volume radius of K is the radius of a Euclidean ball whose 
volume is equal to that of K. Our approach will determine the volume radius 
of E up to a factor which is a power of logci, in particular precisely on the 
scale of powers of d; this is the principal result of this note. The corresponding 
problem for A, the unit ball in the trace class norm, is much better understood. 
Indeed, two-sided estimates for the volume radius of A involving a rather large 
(but universal, i.e., independent of N) constant follow from an early paper |23| . 
Moreover, explicit formulae for the volume of T> involving multiple integrals can 
be produced; see jTH] for an analysis of a closely related problem, which can be 
routinely modified to yield similar expressions for T>. After a preliminary version 
of the present note has been circulated, the author has learned that this circle of 
ideas has led to a closed formula for the volume of T> in a very recent work [28j; 
see Appendix E for more details and |22] f° r related results concerning the Bures 
volume. [Undoubtedly, formulae for the volume of A may be similarly obtained.] 
The unified argument for estimating the volume radii that is presented in this 
paper allows to deduce (from known facts) the value of the volume radius of T> 
up to a factor of 2. 

For the first point, i.e., comparing the volumes of convex sets and their sym- 
metrizations, we use a 1958 result of Rogers and Shephard (see Appendix C for 
more details and background) to deduce that 

2 2 2 n 
-= vol!? < volA < -= vol£>, (4) 

Vd ~ ~ Vd 71 + 1 ' 

where n = dimD = d 2 — 1. [The factor 2 /yd appears because it is the distance 
between the hyperplanes containing T> and — T>; note that strictly speaking we 
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should be writing vo\ n V and vol n+ i A to refer to n and n + 1- dimensional volume 
respectively.] Similarly 

2 2 2 n 
-= vol 5 < volS < -= vol S. (5) 

Vd ~ ~ Vd n + 1 

Combining (J3J) and © we obtain 

2 n V 1 volS vol S 2 n volS 



< < • (6) 

n + 1 J vol A _ vol V ~ n + 1 vol A v ' 

Given that the proper homogeneity is achieved by raising the volume ratios to 
the power 1/n (or l/(n + 1)), we see that one may replace T> and S in (J2J) by A 
and S with the accuracy of the estimates affected at most by a factor of 2. 

It remains to estimate vol A and volS; this will be accomplished by compar- 
ing each of these bodies with the (i 2 -dimensional Euclidean ball Bhs (the unit 
ball with respect to the Hilbert-Schmidt norm; we shall also denote by Shs the 
corresponding d 2 — 1-dimensional sphere). 

Concerning A, we claim that its volume radius satisfies 

1/Vd < (volA/vo\B HS ) 1/d2 < 2/Vd (7) 

To show this, we note first the "trivial" inclusions Bns/Vd C A C Bhs, which 
just reflect the inequalities || • H2 < || • ||i < Vd\\ ■ H2 between the trace class and 
the Hilbert-Schmidt norms. The first inclusion implies the lower estimate on the 
volume radius in (JJJ). The upper bound is less obvious, but it may be shown 
by the following rather general argument. The first step is the classical Urysohn 
inequality, which in our context asserts that 

i/d : 



( volA V 7 f „ ,„ 



Shs 

where the integration is performed with respect to the normalized Lebesgue mea- 
sure on the Hilbert-Schmidt sphere. (For clarity and to indicate flexibility of the 
approach we shall present a general statement and a short proof in Appendix A.) 
The quantity fi d is most easily handled by passing to an integral with respect to 
the standard Gaussian measure, which reduces the problem to finding expected 
value of the norm of the random Gaussian matrix G = G(u) G AiX, usually called 
the Gaussian Unitary Ensemble or GUE. It is well known that E||G|| op = 7^2/^, 
where 7^ := v / 2r(^-)/r(|) for k E N (this uses just 1-homogeneity of the 
norm), and it is easy to check that \Jk — 1 < 7^ < yk for all k. In other words, 
fid ~ E || G || p/rf for large d. On the other hand, it is a well-known strengthening 
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of Wigner's semicircle law that E||G|| op /-\/rf — > 2 as d — > oo. This shows the 
second inequality in (J7|) with 2 replaced by 2 + o(l). We sketch the argument 
that gives the exact number 2 in Appendix F (it follows from known facts, but 
appears to have been overlooked in the random matrix literature), yet we will not 
dwell on it as it intervenes only in the lower estimate in (j2J) and, in any case, the 
constants in our final results are not meant to be optimal. Indeed, (JJJ) combined 
with (J3J) implies that the volume radius of V is between \d~ l l 2 and 2d~ 1 ^ 2 , while 
the formulae from [2E] allow to deduce that it is equivalent to -Jjid^ 1 ^ 2 as d — > oo. 

We now pass to the analysis of the the volume radius of E. We shall show 
that 

l/d 1+a < (volS/volfi^5) 1/d2 < C^\ogd\og\ogd/d 1+a , (9) 

where a is the same as in (J2J). Our main result (J2J) follows then by combining (JTj), 
(Q and ©. [To be precise, one obtains a priori 1/d 2 in the exponent, but a more 
careful analysis of lower order factors such as 2 /yd and l/(n + 1) appearing in 
(jlj)-© allows to replace d 2 by diniiS = d 2 — 1 without any loss in the constants.] 
Before proceeding, let us compare (0) with the results of [6 , which estimate 
from below the in-radius of S in the Hilbert- Schmidt metric by a quantity that 
is of order of d~ v , where r] = 3/2. The easy upper bound on that radius is 
the in-radius of T>, which equals 1/ ^d(d — 1) = 0(d^ 1 ). The second inequality 
in JUJ) yields (for large N) a better upper estimate that roughly corresponds to 
rj = 1 + a ~ 1.094361; we elaborate on these and related issues in Appendix H. 

The first step towards showing will be to replace the sets T> N by their affine 
images which are more "balanced;" this will also explain the appearance of the 
mysterious number a in the exponents. 

Consider first the sets in question when N = 1. As is well known, Si and T>\ 
both coincide with the Bloch "ball," which geometrically is a (solid) Euclidean ball 
of radius 1 / a/2, the boundary of which is the Bloch sphere T\ consisting of pure 
states on £>(C 2 ) (further identifiable with rank 1 projections on C 2 ). Accordingly, 
Si = Ai is a 4-dimensional cylinder whose base is the Bloch ball and whose axis is 
the segment [— J 2 /2, I2/2] of Euclidean length a/2. For definiteness, let us identify 
Aif 1 with M 4 via the usual basis {/ 2 /a/2, a x j a/2, a y j a/2, a z j a/2}, where a x , a y and 
a z are the Pauli matrices (the factors 1 / y/2 make this basis orthonormal in the 
Hilbert-Schmidt sense). Let now A be a linear map on Adf? which is diagonal 
in that basis and whose action is defined by AI 2 = I%/y/2, Aa^ = v/ 3 / 2 cr, for 
i — x,y, z. Set E x := AEj; the important properties of A and Ei are 

(i) the image of the Bloch sphere ATi =: Ti is geometrically a 2-dimensional 
sphere of radius v^3/2 centered at I2/V8 and, as the Bloch sphere itself, it is 
contained in the unit Euclidean sphere Shs] this implies that Ei C Bus 
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(ii) det A = ^27/16 and so vol Si = ^27/16 vol S a 

(iii) vertices of any regular tetrahedron inscribed in 7{ form an orthonormal basis 
in M S 2- 

The geometric property of the set Si, which arguably is the reason for its rele- 
vance, is that the ellipsoid of smallest volume containing it (the so-called Lowner 
ellipsoid of Si) is the Euclidean ball. An equivalent and perhaps more natural 
point of view would be to compare Si with its own Lowner ellipsoid. This is in 
turn equivalent to replacing the Hilbert-Schmidt inner product (u, v) = truv with 

tr((Au){Av)) = (3 truv - tratro)/2. (10) 

It is likely that this non-isotropic inner product and objects associated with it 
play an important role in the theory. In particular, we obtain this way ellipsoids 
which - from the volumetric point of view - are nearly non-distinguishable from 
S or S, and which still enjoy certain permanence relations with respect to the 
action of the unitary group. 

If N > 1, we set S = t N := A 9N V N . Since det A® N = (det A)*" 4 " -1 = 
((27/16) JV / 8 ) d2 = (2 aN ) d2 = {d a ) d \ we deduce that © is equivalent to 

1/d < (vol t/vo\B HS ) l/d2 < CyiogcHoglogd/d, (11) 

For the lower estimate in (fTT|) we shall produce a simple (and seemingly not 
very optimal) geometric argument. Let U\, u 2 , u 3 , tt 4 be vertices of any regular 
tetrahedron inscribed in 7[. By the property (iii) above, (uj)^ =l is an orthonormal 
basis of Aif 1 . Accordingly, the set U := {uj 1 (g) Uj 2 <8> . . . <g> Uj N }, where each ji 
ranges over {1,2,3,4} is an orthonormal basis of Ai s d a - The first inequality in 
ffTTj) follows now from the inclusions U,—U C Sjv and Bns/d C conv(— U U U) 
(the latter is a consequence of the orthogonality of elements of U). 

The above argument may appear rather ad hoc, and so it may be instruc- 
tive to rephrase it in the language of geometry of Banach spaces. Let A\ be 
a linear map from M 4 to M.^ which sends the standard unit vector basis onto 
vertices of any regular tetrahedron inscribed in By construction, A\ is a 
contraction from i\ to C\ and so its iV-th tensor power Af induces a con- 
traction between the respective projective tensor powers of if and C\ (where l\ 
denotes R k endowed with the norm = ^2 As the projective tensor 

product of £i-spaces is again an £i-space, it follows that S^r contains the im- 
age under Af N of the unit ball of if, and hence the image of the Euclidean 
ball of radius 1/d. In particular, vol T l /vo\(Af N B HS ) > (1/d)^. On the other 
hand, one verifies (directly, or by noticing that A = {A^ 1 ] = (A^ 1 * A^ 1 ) 1 ^ 2 ) that 
vo\(Af N B HS ) = ((16/27) N / 8 ) d \o\B HS , which substituted into the preceding es- 
timate gives exactly the first inequality in (jllj) . 
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We note that using for the above calculations the larger volume of the image 
of the if ball (resp., conv(— UUU)) would only result in a slightly better constant 
c in (0) (the value c = ^e/Rir that was mentioned earlier). This is because the 
volume radius of the unit ball in £™ is roughly the same as that of the inscribed 
Euclidean ball, the ratio between the two is ^Jlefix (1 — 0(l/m)). This property 
is behind many striking phenomena discovered in the asymptotic theory of finite 
dimensional normed spaces, and is closely related to our upper estimates for vol S 
and vol S, to which we pass now. 

To prove the upper estimate in (jlljl . we shall again use the Urysohn inequality. 
Analogously to (jHJ) and to the reasoning that followed it, we get 

- \ i/d 2 

— < / max tx{XA) dA = V^E max tr(XG) (12) 

VO\B HS J J Shs x ^ x ^ 

and so it remains to show that the above expectation is O ( \/log d log log d) . The 
expression under the expectation can be thought of as a maximum of a Gaussian 
process indexed by £ (this just means the family tr(XG(u)), X 6 S, of jointly 
Gaussian random variables). There are several methods of differing sophistication 
which can be used to estimate the expectation of such a maximum. The two 
leading ones are the Fernique-Talagrand majorizing measure theorem, which gives 
the correct asymptotic order, but is usually difficult to apply, and the Dudley 
majoration (by the metric entropy integral), which is almost as precise and usually 
easier to handle; see |T2J for a comprehensive exposition. We shall employ here 
an even simpler "one-level-discretization" method which, in our context, yields 
approximately the same result as the Dudley majoration, and which we now 
describe in elementary language. 

Let /i be the standard Gaussian measure on R m (i.e., the one given by the 
density m(x) = (27i)~ m ^ 2 exp(— \x\ 2 /2), where | • | is the corresponding Euclidean 
norm) and let F C M. m be a finite set contained in a ball of radius R. Then 

/ max (y, x) dfi(x) < R^2 log(#F), (13) 
J R m y£F 

where ^ stands for the cardinality of a set. The estimate above is usually quoted 
with a different numerical constant appearing in place of 2, but it is not difficult 
- even if somewhat tedious - to verify that it holds in the form stated above. 
The idea now is to construct a finite set F C X such that convF D r£ for an 
appropriate r G (0, 1); it will then follow that 

Emax tr(XG) < r~ l yj2\og(#F). (14) 
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[We note that the maxima of the type appearing in (|T2"j) . (fT^j) or (JT3J) do not change 
if we replace the underlying (closed) set F by its convex hull or, conversely, by 
its extreme points.] Specifically, F will be a "sufficiently dense" subset of the set 
of extreme points of E, i.e., of —T U T, where 

f = f N = {Ap± <g> Ap 2 <g> . . . <g> Apiv} 

and where each p,- is a pure state on £>(C 2 ) (i.e., an element of the Bloch sphere). 
In other words, T is a tensor product of N copies of F\ = BT X which, as we noted 
earlier, is geometrically a 2-dimensional sphere of radius a/3/2 contained in the 
unit sphere of the 4-dimensional Euclidean space. 

We start by constructing an appropriate dense subset (usually called net) in 
each copy of T\ and then consider tensor products of those nets. To facilitate 
references to existing literature we first look at the unit Euclidean ball S 2 rather 
than 7j. Let S E (0, a/2) and let M be a 5-net of S 2 , i.e., a subset such that 
the union of balls of radius 5 centered at points of M covers S 2 . An elementary 
argument shows that convA/" contains then a ball of radius (1 — S 2 /2) centered at 
the origin. If now F\ C / 7 1 is an appropriate dilation of M (i.e., with the ratio 
a/3/2), then convFi contains a ball of radius (1— 5 2 /2)-a/3/2 (in the 3-dimensional 
affine space containing 71) with the same center as that of of 7[. It follows that 
conv(-Fx U F x ) D (1 — 5 2 /2) E a and consequently if we set F := (-Fx U F ± ) m = 
F® N U then convF D (1 - 5 2 /2) N t N . 

In remains to find a reasonable bound on #F = 2(#F i ) JV = 2(#A/') JV . A 
standard argument comparing areas of caps and that of the entire sphere shows 
that one may have a 5-net of S 2 of cardinality < 16/5 2 . [This bound is not optimal; 
the asymptotically - as 5 — > - correct order for cardinalities of efficient £-nets of 
S 2 is (2/v / 3) 3 vr/<5 2 , but coverings of Euclidean spheres, even in dimension 2, do not 
appear to be completely understood.] This leads to an estimate #F < 2(16/<5 2 ) 7V , 
which in combination with (pPf]) gives 



Emax tr(AG) < (1 - 5 2 /2)~ N a/2 log(2(16/<5 2 H. (15) 

Optimizing the expression on the right hand side over 5 G (0, a/2) yields a quantity 
that is of order y/2N log iV for large iV (choose, for example, 5 = (NlogN)^ 1 ^ 2 ), 
as required to complete the proof of © (and hence of ©)• Moreover, substituting 
the obtained bound into ()12j) and verifying numerically small values of N yields 

i/d 2 



volE \ < y/4iV log 2 (4AQ _ A/41og 2 tZ log 2 (4 log 2 jjj 



voIBhs I d d 
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(note that the inequality is trivial for N = 2), which implies that (J2J) holds with 

the third member of the form ^ ^i/° 2 + 2 i — ^ ■ This may be somewhat improved for 
small to moderate values of N by using estimates on cardinalities of nets of S 2 
listed in |3J. see Appendix G. 
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Appendix A The Urysohn inequality. If K is a convex body in the Tri- 
dimensional Euclidean space which contains in its interior, then 

where \\x\\k is the gauge of K (the norm for which K is the unit ball if K is 0- 
symmetric - which is the case in the main text) and the integration is performed 
with respect to the normalized Lebesgue measure on the sphere S m ~ l . If K is 
0-symmetric, this may be combined with the Santalo inequality [20] which asserts 
that 

vo\K vo\K° < 
voW ' vo\B ~ ' 

where K° := {x : (x,y) < 1 for all y G K} is the polar of K, to obtain 
(volA^/voLB) 1 /" 1 < L , 1 1 a; 1 1 ^ cfe, and the Urysohn's inequality 

(vo,A'/vo. B )V™< f M *-d*=f ^( X ,y) dx 
Js™- 1 Js™- 1 yeK 

follows by exchanging the roles of K and K°. Moreover, since the Santalo inequal- 
ity holds also for not-necessarily-symmetric sets after an appropriate translation, 
the above inequality holds for any (measurable) bounded set K. Indeed, the in- 
tegral on the right equals 1/2 of the mean width of K, a well-known classical 
geometric parameter of a set in the Euclidean space, which does not change if K 
is replaced by its translation. It is primarily the mean width of various sets - and 
not directly volume - that is being majorized throughout this paper. 



10 



This is not the most elementary proof of the Urysohn inequality, but one that 
offers a lot of flexibility. For example, the repeated applications of the Holder 
inequality in the first chain of inequalities above can be modified to yield as the 

last expression (J S m-i \\ x \\ P k^ x ) ^ P f° r an arbitrary p > and, letting p — > 0, 
the geometric mean exp f— J Sm _i log ||;r||i<:(ia;). Similar inequalities also hold if 
p G [—n,0), and the case p = —n is of course the strongest statement of such 
nature, the Santalo inequality itself. We also take this opportunity to point 
out that the Santalo inequality and the so-called reverse Santalo inequality j2] 
together imply that the volume radius of a convex set and its polar are roughly 
(i.e., up to universal multiplicative constants) reciprocal. 

The application of the Urysohn inequality in (JHJ) uses implicitly the elementary 
fact that || • ||a° = || • || p- in other words, the trace class norm || • ||i and the 
operator norm || • || op are dual with respect to the trace duality. For the sets A at, the 
Urysohn inequality gives the correct order of the volume radius, cf. ((7j). However, 
this is not always the case, even for rather regular convex bodies. For example, if 
K is the unit ball of i™, then its volume radius - as we have already mentioned 
- exceeds the radius of the inscribed Euclidean ball by less than *J2epK, while 
the upper bound obtained from the Urysohn inequality contains a parasitic factor 
which is of order \J\ogm. It is thus conceivable that the logarithmic factors in 
can be replaced by universal numerical constants. On the other hand, if we 
do use the Urysohn inequality to establish an upper bound for the volume radius 
of £ (cf. ffTTj) . (fT2j) . (JHJ)), then our estimates can not be substantially improved. 
Indeed, since our argument showed that £ contained a rotation of the unit ball of 
if , the previous remark implies that the = ^J\og^d factor can not be then 
avoided, and so it is only the -y/log log d factor that can possibly be eliminated by 
more careful majorizing of E||G||£o. 

Finally, we mention that there exist general volume estimates for convex hulls 
of finite sets which are asymptotically more precise than the one we derive from 
the Urysohn inequality (see pQ and its references). However, these estimates 
are equivalent to the ones presented here in the relevant range of parameters 
and, moreover, their formulations available in the literature contain unspecified 
numerical constants, which would make the corresponding bounds difficult to 
apply for specific values of N . 

Appendix B Projective tensor products of normed spaces. If X and Y are (say, 
finite dimensional) normed spaces, their tensor product X £g> Y may be endowed 
with the projective tensor product norm || • 1 1 ^ defined by 

IMU := inf {E7=i INI ' 11%' II : EJU x i ® Vi = T } ■ 
The resulting normed space is usually denoted X®Y or X ® n Y . If Bx and 
By are unit balls of X and Y respectively, it follows that the unit ball of X®Y 
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coincides with 

conv {x ® y : x G B x , y G B Y } = conv {x <g> y : x G exti?x, 1/ G ext-By }, 
where ext(iT) denotes the set of extreme points of K. If X = Y = Cf, the 
analysis is further simplified by the fact that the set of extreme points of A 1; 
the unit ball of Cf, is of the form —Ti U where T\ is the set of pure states. 
The fact that the set S 2 is the unit ball of C\®C\ follows directly from these 
identifications. Projective tensor products of more than two spaces are defined 
analogously (or by induction), and one similarly checks that the unit ball of the 
iVth projective tensor power of Cf is Sjy. Tensor products involving spaces C k 
with k > 2 may be treated in the same way. For example, the symmetrization 
of the set of separable states on B(C kl (g> C k2 <8> • • .C hm ) can be identified with 
the unit ball in C\ 1 ®C X 2 <8> ■ ■ ■ ®C k m . The problem of the relative size of the set of 
separable states on B((C D )® N , or iV qudits, was investigated in [T^]. While a more 
definitive treatment of the higher- dimensional case will be presented elsewhere, 
we offer some comments on the topic in Appendix I. 

On the more elementary level, the projective tensor square of a Euclidean 
space C k <g)C k can be identified with (A4k, \\ ■ ||i) and the contractively comple- 
mented subspace of its self-adjoint elements is, in our notation, C\. We also 
mention in passing that the dual space to X®Y can be identified with the so- 
called injective tensor product of the duals X* and Y*, and so S° can be thought 
of as a unit ball in the injective tensor power of the Hermitian part of B(C 2 ). 
While making this identification explicit doesn't seem to help our analysis at the 
present level of depth, we mention in passing that various existing criteria for 
detecting entanglement use separation theorems for convex sets which are based 
on a form of this duality. 

Appendix C The Rogers- Shephard results on symmetrizations of convex sets. 
Let W C M n+1 be an n-dimensional convex set and denote by h the distance from 
the affine hyperplane H spanned by W to the origin. Let Q be the symmetrization 
of W, i.e., Q := conv(— W U W). It was shown in ^7] that then 

2h volW < volQ < 2h vo\W, (17) 

n + 1 y ' 

where by volVT and volfi we mean the n- and the n + 1-dimensional volume 
respectively. To explain the factors appearing in (fTTj) we note that the inequalities 
becomes equalities if W is centrally symmetric for the first one (this is simple) 
and if W is a simplex for the second (this is the heart of the Rogers-Shephard 
result). 

To further clarify the first inequality in (|17j) (used in the upper estimates on 
the volume of separable states, which is the main point of this note) we point 
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out that it is actually a simple consequence of a much older theorem of Brunn- 
Minkowski, and more specifically of the following corollary of that theorem. 

Let K be an n + 1- dimensional convex body, u - a vector in the ambient 
space containing K and H - a hyperplane in that space. Then the function 
t — > vol(i^ fl (tu + H)) l l n (the n-dimensional volume) is concave on its support. 

If we apply the above fact with K = Q and u - a unit vector perpendicular 
to H, then the function <p(t) := vol(fi fl (tu + H)), being even on [—h, h], must 
attain its maximum at and minimum at h and —h. The first inequality in 
(fTTjl follows then from the Cavalieri principle. A version of the second inequality, 
which would be sufficient for our purposes, follows similarly from the estimate 
(f)(0) < 2~ n ( 2 ™)vo\W, which is the main result of 

Appendix D Working directly with non-symmetric sets. Similar but slightly 
more complicated arguments may be used to obtain upper estimates for the vol- 
umes of the non-symmetric sets T> and S by studying directly these sets and not 
their symmetrizations £ and A. [In principle, this could help to avoid the para- 
sitic factors 2 n / (n+ 1) - where n = d 2 — 1 - when passing from S, A to S, T>.] In 
both cases it is convenient to pass to a translate of the set in question obtained by 
subtracting the appropriate multiple of the maximally mixed state Id/d, and to 
consider the translates as subsets of H , the d 2 — 1-dimensional space of matrices 
with vanishing trace. 

For the set V (translated by Id/d), the quantity which replaces || • || op in the 
analogue of (JHJ) is Ai(-), the largest eigenvalue of a matrix. This is of course 
dominated by the norm, and since the (random Gaussian) trace matrix Go 
can be represented as a conditional expectation of the general Gaussian matrix 
G, it follows - by the convexity of the norm or of the largest eigenvalue - that 
EAi(Go) < M G \\op < 2v^ which, after some work, leads to an upper estimate 
for the volume radius of T> identical to that of A obtained in (J2J). [To fully justify 
the steps above one needs to appeal to Appendices A and F.] 

For the set S, we pass first to the face S of the rescaled set £ that corresponds 
to S, and then subtract Id/d 3 / 2 (the difference with respect to the case of T> is 
due to the rescaling). Next, we "approximate" the translate by sets built from 
the points corresponding to elements of F. There are several differences between 
this setting and that described in the main text, but they can be accounted for 
fairly easily. The good news is that the new points are not on the sphere since the 
component in the direction of Id was subtracted, but this improves our estimate 
on vol S only by a factor 1 — 0(d~ 1 ). A loss which is even less significant is due to 
the fact that by reducing the dimension by 1 our formulae will involve the quantity 
7rf2_! rather than 7^2. A somewhat more substantial loss comes from the fact that 
due to the rescaling the width of £ in the direction of Id is different from that 
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of E by a factor of 2 N I 2 = d x l 2 \ this affects the relationships between volumes of 
these bodies and those of S and S, and consequently our estimates, by the same 
factor. However, since we are in dimension d 2 — 1, the loss in the volume radius 
is a not-so-significant 1 + O (log d/d). The final issue that needs to be analyzed 
is that while we knew that convF D (1 — 5 2 /2) N H, it is not a priori clear 
that a similar inclusion holds for the face S (or, more precisely, for its translate 
iS — Id/d 3 / 2 ). While for a general convex set K C H\ the relationship between K 
and its symmetrization conv(— K U K) may be more involved, using the fact that 
S is a convex hull of points contained in a sphere centered at Id/d 3 ^ 2 we can infer 
that convF D (1-k)E implies conv(FnS) D / d /d 3/2 + (1 - 2k) (5 - I d /d?l 2 ) } and 
the difference between the factors 1 — k and 1 — 2k is asymptotically insignificant. 

While the above argument allows to avoid the symmetrizations while estimat- 
ing the volume radii of T> and S from above, we get the same majorants as in (J7J), 
0. Moreover, we still need to look at the symmetrizations for the lower bounds 
which are needed to derive (J2J) . [It is also possible to obtain a lower bound on the 
volume radius of S directly by noticing that S contains a simplex spanned by the 
elements of (A® N )~ 1 U, but since for the simplex the relevant Rogers-Shephard 
inequalities become equalities, there is again no significant improvement.] 

Appendix E The exact expressions on the volume of T> and the constants in 
(J2J) as N — > oo. After a preliminary version of this note has been posted, the 
author has learned that a closed formula for the volume of T> was derived in a 
very recent work [2Bj- While we were able to calculate the volume radius of T> to 
within a factor of 2 by the same methods that were employed to analyze S and 
with very little extra work, it is instructive to compare the so obtained estimates 
to those that can be deduced from the exact formula, which in our notation reads 

volP(C d ) = Vd (2tt)^- 1 )/ 2 (18) 

T{d 2 ) 

A tedious but routine calculation based on the Stirling formula shows that the 
volume radius of V(C d ) behaves as -Jjjd~ l / 2 {1 + 0(d~ 1 )) as d — > oo 

We now recall the refinements related to the volume radii of Sn and 
suggested in the main text. First, we had the argument that gave c = ^/e/8ir 
in (J2J), based on using the volume of the unit ball in if, i.e., 2 d2 /(d 2 )\, as a 
lower bound for volE (see the comments following (|ll|l). Next, we noted that, 
for large N, the expressions in ()15|) can be majorized by a quantity that is of 
order y/2NlogN. Combining these with the improvement related to T>^ we 
are led to an asymptotic version of © with c N -> e 3 / 4 /V27r « 0.844561 and 
C N -> e 1/4 v/2/log2 w 2.1811. 

Appendix F Norms of GUE matrices and the constant 2 in (|7j). It has been 
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known for some time (in fact in a much more general setting) that if G = G(u) 
is the random matrix distributed according to the standard Gaussian measure 
on Ai™ (usually called the Gaussian Unitary Ensemble or GUE), then, for large 
d, \\G\\ op is, with high probability, close to 2\/d . We sketch here a derivation, 
from known facts, of the arguably elegant inequality E||G|| op < 2yd, valid for 
any d, which appears to have been overlooked in the random matrix theory liter- 
ature. Similar inequalities are known for Gaussian matrices all whose entries are 
independent or for real symmetric matrices (also known as the GOE ensemble; 
however, in the latter case the precise inequality seems to have been established 
only for the largest eigenvalue, and not for the norm), see [3]. Analogous inequal- 
ities with the expected value replaced by the median can probably be deduced - 
at least for large d - from fl%\ EE] . 

Our starting point are the recurrence formulae for the (even) moments a p = 
Op(d) := rf" 1 Etr((G/2) 2p ), pGN, derived, e.g., in [7] (see also |HJ, formulae (6) 
through (9), for a similar argument and a related estimate) 



2p 



I / p(p - 1) 2p - 3 

a P-! H 732 nZ a P-2 



2p + 2 \ " 1 ' Ad 2 2p 
with ao = 1 and a± = 1/4. From these one easily derives by induction 



1 



ap ~ 2* P (p + l) 



2p 
P 



n(i+ 

3=1 



j(j - 1) 

Ad 2 



[This estimate is actually asymptotically precise for p = o(d).} Next, using suc- 
cessively the Stirling formula to majorize the binomial coefficient, the inequalities 
1 + x < e x and Y7j=iiU ~ 1) — P 3 /3 to estimate the product, and denoting 
t = pd~ 2 ^ 3 , we arrive at 



Etr((G/2) 2p ) =da p <d 



3 p 3 /12d 2 



Tip 



3/2 



e t3 ' 12 

y^tfl 2 ' 



Hence 



l -nG\W< (Etr((G/2) 2 ^)) 1/2p < 



l/At 



1/d 



2/3 



This is valid for t > 0, at least if the corresponding value of p = td 2 ^ 3 is an 
integer. The minimal value of the expression in brackets over t > is attained at 
t w 1.38319 and is approximately 0.738542 w exp(-0.303077) < e" 03 . Since for 
sufficiently large d the interval corresponding to values which are < e~ ' 3 contains 



15 



an element of <i _2 / 3 N, we deduce that for such d we have E||G|| op < 2e _a3d 2/3 . 
A more careful checking shows that in fact the inequality E||G|| op < 2 — 0.6<i _2//3 
holds for all values of d (in fact, by the above argument, the same upper estimate 
is valid for (E||G||£ J 1 / 7, with, say, r = 2 or r = d 2 ^). 

Going back to the issue of having the precise constant 2 in inequality (JTj), let 
us note that the other source of difficulty, namely the fact that the parameter 7^ 
is only asymptotically of order \fk but not equal to y/k, introduces an error that 
is of smaller order than our "margin of safety." As pointed out earlier, we have 
7fe > \Jk — 1 and so 'jk/Vk > \Jl — 1/k) ~ 1 — l/2k. The relevant value of k is 
d 2 , leading to the relative error of order d~ 2 /2, as opposed to the margin of safety 
of 0.3d~ 2 / 3 (note also that 4 is the smallest value of d that is of interest). 

Appendix G Upper estimates on voISn for small to moderate N. We now 
indicate how one may use the explicit efficient nets of the sphere S 2 listed in jS] 
to majorize the volume radius of Sn if N is not too large. As a demonstration, 
we will derive bounds for the volume of the set of separable states on 8 qubits 
(one may say, a qubyte). 

The site jH] lists, for m € {4, ... , 130}, sets M m of m points in S 2 such that 
every point of S 2 is within e = e m (measured in degrees) from one of the points of 
Am, with the dependence m — ► e m "putatively optimal" (and very likely nearly 
optimal). Noting that 5 = 2 sine/2 we verify numerically that in most of the 
interesting range the putatively optimal value 5 m verifies ~ 5 (more pre- 

cisely, 5.1 ± 1%, still not far from the asymptotic value (2/y / 3) 3 7r w 4.8368 that 
we mentioned earlier). Since in the present context the bound in (|15jl becomes 
(1 - 5 2 m /2)- s ^j2 log(2m 8 ) = (cose m )"V 2 log(2m 8 ), substituting m = 5/5 2 a leads 
to a function (f)(5) = (1 — 5 2 /2)~ 8 a/2 log(2(5/<5 2 ) 8 ), which attains its minimum 
very near 5 = .15, which suggests that the optimal value of m should be around 
222. This is beyond the range of the tables from [H], but using the largest avail- 
able m = 130 and the corresponding €130 = 11.3165625° yields an upper bound of 
10.417406, which is less than 2% larger that the majorant that would presumably 
be given by m = 222. Plugging in the obtained bound into (j!2)l and using the 
fact that (volS/volS) 1 ^ 2 = d a = (27/16) 7V/8 we are led to 

v := (vol£/voLBtf 5 ) 1/d2 < (16/27) • 10.417406/7^ < 0.02411446 

Taking into account (jSJ) and substituting the explicit expression for voWhs we 
obtain 

,„ Vd n d2 / 2 v d2 

vols - -T rK/2 + i) ' 
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Finally, using the closed formula for voW from (|18|) we get 

volS < (2n)^ v d2 T{<?) < g g 1Q _ 19996 

void - 2 d2 /2 +1 r(d 2 /2 + 1) njLi r(j) 

which (modulo rounding errors) is equivalent to a much less impressive bound 
(volS/vol£>) 1/dimI? < 0.49534. 

The smallest N for which an argument such as above gives a non-trivial upper 
bound appears to be 6; we get (vol56/ / vol£ , 6) 1//dim ' DB < 0.95. We refer to [21] and 
its references for extensive (largely numerical) treatment of the case N = 2. 

Appendix H The in-radii oJSn and Ejv- The papers [HE] estimate from below 
the in-radius of S in the Hilbert-Schmidt metric by a quantity that is of order of 
d~ v , where rj = log 20/ log4 « 2.160964 and 3/2 respectively. The "trivial" upper 
bound on that radius is the in-radius of T>, which - by a rather elementary and 
well-known argument - equals 1/a/ d(d — 1) = 0(d^ 1 ). By comparing volumes 
we see that the second inequality in (j^J) yields an asymptotically better upper 
estimate that (up to logarithmic factors) corresponds to rj = 1 + a ~ 1.094361. 
This follows by taking into account (jSJ) or by observing that, by simple geometric 
considerations, the Euclidean in-radius of £ is at least as large as that of S (the 
latter considered in the hyperplane H\ of trace one matrices). By tinkering with 
the argument it is possible to remove the logarithmic factors and, indeed, to 
obtain an upper bound on the in-radius of £ (and hence of S) which is o(<i -1 ~ a ), 
but to improve the exponent new ideas would be necessary - if that is at all 
possible, that is. 

Let us also note that our argument yields as well a lower bound 6~ N / 2 on the 
in-radius of £ which corresponds to rj = log6/log4 « 1.3863, and so is stronger 
that those that can be formally derived from [3] or To see this it is enough 
to combine the "trivial" lower estimate d~ l = 2~ N on the in-radius of the set Si- 
delined in what follows (cf. the paragraph following (jll|) ) with the known value 
(3/2) N / 2 of the norm of the related map A® N . However, since the in-radius of 
S may a priori be (and probably is, except when N = 1) smaller than that of 
S, this doesn't improve the lower bounds from jHIH]- It is also likely that the 
methods from those papers may yield a lower bound on the in-radius of £ that 
is better than 6~ N / 2 . 

While our calculations narrow down the potential range of the in-radii of S 
and E, and while further progress along the same lines may be possible, it seems 
likely that to obtain the exact values of exponents a more careful calculation 
involving, e.g., spherical harmonics may be necessary. On the other hand, as we 
have already noted, it may be more natural to consider in this context the inner 
product norm, which is different from the Hilbert-Schmidt norm and induced by 
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the inner product (u, v) — *> (3trm> — trutru)/2 on each factor Aif 1 . The (solid) 
ellipsoid S, which is the unit ball with respect to the corresponding inner product 
norm on Ai s d a = , verifies S/d C S C S (this is just a restatement 

of Busjd C S C Bus) an d these inclusions are essentially optimal. On the 
one hand, every pure state clearly belongs to the boundary of 8. On the other 
hand, it follows from, say, (fTTj) that the volume radius of S is, on the power 
scale, the same than that of S/d. This means that - from the volumetric point 
of view - E and S/d are nearly non-distinguishable. It would be of interest 
to determine the precise in-radius of S with respect to the inner product norm 
induced by S or, equivalently, the largest e such that Id/d + e(SnH ) C S, where 
H := {A G M. s d a : tiA = 0}. [It is conceivable that that radius is of order d -1 .] 

Appendix I Separable states on N audits. (See also t 19j.) Most of the elements 
of our analysis can be generalized to tensor products involving spaces C D with 
D > 2, leading to non-trivial but not definitive results. As a demonstration, 
let D > 3 and consider the family of spaces Ti^ = (C D )® N . We shall employ 
analogous notation to that of the main text, in particular d = dimTi,N = D N . 
The set of pure states on £>(C D ) coincides with the projective space CP ' 1 whose 
real dimension is 2D — 2 and which admits, for 5 > 0, 5-nets of cardinality not 
exceeding (C'/5) 2D ~ 2 , where C is a universal constant. This leads to a bound 
on (voIEjv/voLB^s) 1 ^ 2 which is of order (1 - 5)~ N y/ND log {C'/S)/d. Choosing, 
say, 5 = 1/N and using the same bound on voWn as earlier combined with the 
"easy" part of (0) we obtain 

(vol^/volPjv) 1 /^ = O (y/ND logiV/d 1/2 ) . 

Since d = D N , this leads to a non-trivial bound even for iV = 2 if D is large 
enough. It is also possible to improve slightly the exponent of d by working (as 
we did for D = 2) with a more balanced affine image of Sjv- The resulting 
improvement « = «o decreases as D increases; for example, = f-j^|f — \ ~ 
0.0608264 and, for large D, old ~ {2D log D)~ l . However, showing optimality of 
the so obtained exponents requires - for D > 2 - new ideas and will be presented 
elsewhere. 
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